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ABSTRACT 


The  S-matrix  of  a  many-body  system  is  usually 
assumed  to  be  an  analytic  function  of  the  coupling 
constant.  A  system  of  quasi-particles  interacting 
through  a  pairing  force  is  considered  .  It  is  shown 
that  for  this  sytem  the  S-matrix  is  a  non-analytic 
function  of  the  coupling  constant.  This  suggests  that 
the  many-body  problem  should  not  be  treated  using  the 
conventional  perturbation  methods.  Bogoliubov ' s ,  non- 
perturbative ,  operator  expansion  is  generalized.  This 
generalized  quasi-particle  expansion  is  derived  from 
first  principles  and  the  subsidiary  conditions  that  are 
imposed  on  the  expansion  are  discussed.  The  resulting 
lowest  order  equations  are  solved.  It  is  found  that  a 
distribution  function  of  bosons  can  be  defined.  This 
function  depends  on  the  form  of  the  interaction  and  under 
certain  condition  yields  the  Bogoliubov' s  results.  For 
separable  interaction  the  equations  are  solved  exactly 
and  it  is  shown  that  the  energy  spectrum  has  the  correct 
behaviour  in  the  limit  of  very  long  and  very  short  wave 
lengths.  To  investigate  the  life-time  of  quasi-particles 
solvable  model  Hamiltonians,  with  effective  interaction 
between  quasi-particles,  are  proposed.  This  effective 
interaction  can  be  found  from  the  elastic  scattering  of 
two  quasi-particles,  using  the  S-matrix.  It  is  found  that 


. 
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in  the  Born  approximation  the  effective  interaction  is 
density  dependent.  Finally  a  'pairing  type'  of  boson 
operator  expansion  is  proposed.  This  expansion  is 
especially  suitable  for  treating  scattering  processes 
from  a  hard  core  potential.  The  resulting,  G-matrix, 
integral  equation  is  discussed.  It  is  shown  that  the 
effective  interaction  is  just  the  first  order  approxi¬ 
mation  to  the  G-matrix.  A  procedure  for  calculating  the 
energy  levels  and  the  life-time  of  quasi-particles  is 
developed. 
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CHAPTER  1.  INTRODUCTION 


In  1908  Karmerlingh  Onnes  succeeded  in  producing 
for  the  first  time  liquid  helium.  When  some  of  the 
strange  properties  of  the  liquid  became  known,  the  research 
into  the  many-body  problem  started  in  earnest.  As  early  as 
1938  some  theoreticians  began  to  realize  that  liquid  helium 
could  not  be  understood  in  the  same  terms  nor  by  using  the 
same  methods  as  applied  to  ordinary  liquids.  Since  then 
the  theory  of  liquid  helium  has  been  developed  on  three 
different  levels. 

The  first  level  is  the  phenomenological  theory  of 
1  2 

London  and  Tisza  .  Here,  a  comparison  with  the  conden¬ 
sation  of  a  free  Bose-Einstein  gas  leads  to  a  two-fluid 
model  and  predicts  a  macroscopic  occupation  of  a  single 
quantum  state,  the  ground  state,  which  is  identified  with 
the  superfluid  particles.  With  this  analogy  the  unusual 
transport  properties  of  He  II  as  well  as  the  A-transition 
of  liquid  He  are  well  explained. 

The  second  level  is  based  upon  the  semiphenomenolo- 
gical  excitation  picture  derived  from  quantum  hydrodynamics 

3 

by  Landau  .  Here  it  is  assumed  that  the  properties  of 
liquid  helium  may  be  computed  by  treating  it  as  a  gas  of 


. 
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weakly  interacting  elementary  excitations,  the  phonons  and 
the  rotons.  The  energy  spectrum  is  fitted  from  experiment 
and  the  thermal  properties  of  the  liquid  are  calculated 
from  it. 

The  third  level  attempts  to  explain  liquid  helium 

from  first  principles,  i.e.,  on  microscopic  basis.  Here 

however,  we  are  faced  with  the  impossible  task  of  trying 

23 

to  treat  exactly  an  ensemble  of  ~10  helium  atoms  inter¬ 
acting  via  an  appropriate  potential.  In  view  of  this,  in 

seeking  an  adequate  microscopic  theory  to  describe  such  a 

4 

system,  we  should  bear  in  mind  Feynman's  comment  about 
the  limitation  of  such  an  approach:  "The  quantum  mechanics 
will  not  supplant  the  phenomenological  theories.  It  turns 
out  to  support  them."  Thus  in  spite  of  the  great  progress 
made  in  this  direction  by  the  use  of  powerful  and  ingenious 
mathematical  techniques,  the  results  obtained  are,  at  best, 
crude  approximations. 

The  primary  aim  of  every  microscopic  theory,  which 
describes  a  given  system,  is  to  calculate  the  energy  level 
spectrum  of  that  system.  From  this  spectrum  physical  pro¬ 
perties  such  as  pressure,  specific  heat,  compressibility, 
etc...,  can  be  deduced.  Therefore  it  is  essential  that  we 
know  how  to  calculate  this  energy  level  spectrum.  For 
complicated  systems  this  is  indeed  a  difficult  problem. 
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However,  from  experimental  data  we  usually  find  that  in  cer¬ 
tain  energy  regions  the  system  can  neatly  and  conveniently 
be  described  by  elementary  excitations.  The  reason  for 
this  is  that  no  matter  what  the  statistics  of  the  system 
are  and  without  necessarily  assuming  weak  interactions 
between  the  original  (bare)  particles,  the  Hamiltonian 
of  the  system  can  be  reduced  to  the  form 

H  =  Hqp  +  Hqp  ,  (1.1) 

o  I 

where  Hqp  is  the  free  part  of  the  Hamiltonian  describing 
free  (dressed)  particles  (quasi-particles) ,  and  Hqp  des¬ 
cribes  the  interaction  between  the  particles  (quasi-particles) . 
The  reduction  (1.1)  is  by  no  means  unique.  However,  for 
(1.1)  to  be  useful,  Hqp  must  be  small. 


The  advantage  of  the  reduction  (1.1)  is  that  now  the 
properties  of  Hqp  can  be  determined  exactly  while  Hqp  can 
be  treated  as  a  perturbation.  Consider  the  many-body 

5 

Hamiltonian  which  in  the  frame  of  second  quantization 
takes  the  form 


H  =  l  (ek-  U)a 
k  ~ 


kak  + 


1 

j  £mn 


t  t 

6  .  ,  n  ,v  ( j-m)  a  .ana  a 
j+£,m+n  ~  ~  j  i  n  m 


(1.2) 


o 

Here  =  k  /2m  is  the  kinetic  energy,  y  is  the  chemical 
potential,  v(j)  is  the  Fourier  transform  of  a  two  body 
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central  potential  and  a£  and  a^  are  the  creation  and  anni¬ 
hilation  operators  respectively.  We  note  that  the 
Hamiltonian  (1.2)  is  already  in  the  form  (1.1).  Since 
we  know  all  the  properties  of  the  free  part,  Hq ,  of  (1.2) 
we  would  like  to  treat  H^=  H  -  Hq  are  a  perturbation.  Here 
however  one  is  faced  with  diverging  terms  in  the  perturba¬ 
tion  series  which  arise  from  the  macroscopic  occupation  of 
the  zero  momentum  state.  These  so  called  "dangerous 

7 

diagrams",  which  describe  the  self-energy  of  the  system, 
appear  as  a  consequence  of  the  considerable  energy  necessary 

g 

to  excite  particles  from  the  ground  state.  Bogoliubov  in 
/  1947  showed  how  one  can  overcome  this  difficulty.  In  his 
theory,  Bogoliubov  assumed  macroscopic  occupation  of  the 
ground  state  and  proceeded  to  separate  out  the  zero  momentum 
states  from  the  Hamiltonian  thus  eliminating  the  troublesome 
"dangerous  diagrams".  The  resulting  Hamiltonian  was  then 
diagonalized  and  the  energy  spectrum  obtained  agreed  with 
the  observed  results  for  a  weakly  interacting  boson-gas. 

Aside  from  the  semiphenomenological  beautiful  work 
9  10 

of  Feynman  and  Feynman  and  Cohen  ,  most  of  the  work  that 

followed  Bogoliubov' s  theory  was  centered  around  his  main 

assumptions  and  methods.  Penrose  and  Onsager"^  were  first 

to  show  that  Bose-Einstein  condensation  is  possible  in  the 

12 

general  case  of  interacting  bosons.  Beliaev  and  Hugenholtz 
13 

and  Pines  independently  extended  the  Bogoliubov  treatment 
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to  higher  order  approximations  by  applying  the  Green's 

function  formalism  of  quantum  field  theory  to  the  many- 

body  interacting  boson  system.  A  different  approach  was 

14 

proposed  by  Lee,  Yang  and  Huang  who  used  a  hard  sphere 
model  Hamiltonian  with  a  pseudopotential,  which  was  diago¬ 
nalized  by  the  Bogoliubov  transformations.  Brueckner  and 
15 

Sawada  introduced  a  reaction  matrix  or  'T-matrix', 
which  removed  the  difficulties  arising  from  the  singularity 
in  the  interparticle  potential.  Luban's^  pairing  Hamiltonian 

1  7 

and  the  variational  approach  of  Girardeau  and  Arnowitt 

“I  O 

and  Valatine  and  Butler^0  yielded  essentially  the  same 
^results  as  Bogoliubov' s  with  perhaps  better  mathematical 
justification  and  hindsight. 

All  the  theories  we  referred  to  (and  many  more  that 
we  did  not)  have  one  thing  in  common;  it  is  very  difficult 
to  find  or  calculate  higher  order  corrections  to  the  energy 
spectrum  in  a  systematic  way.  While  the  perturbative 
theories  give  us  formal  prescriptions  by  which  we  should 
be  able  to  calculate  higher  order  terms,  these  theories  must 
be  tailored  for  the  treatment  of  superfluids  so  as  to  get 
rid  off  divergences  and  other  difficulties  that  might  occur. 

A  difficulty  of  another  nature  which  is  usually  over¬ 
looked  is  the  question  of  the  convergence  of  the  perturba¬ 
tive  method.  Thus  if  our  interaction  Hamiltonian  H^  is 


. 
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proportional  to  some  parameter  g,  and  the  resulting  power 
series  in  g  that  we  obtain  from  perturbation  theory  con¬ 
verges,  then  we  have  a  good  physical  description  of  the 
system.  Thus  even  for  non-linear  systems  such  as  in 
many-body  theory,  one  always  assumes  that  the  solution 
would  be  an  analytic  function  of  the  coupling  constant, 
for  small  enough  values  of  this  constant  (the  reason  for 
this  assumption  is  discussed  at  the  beginning  of  Chapter 

g 

2).  The  Bogoliubov  theory  of  weakly  interacting  bosons, 

19 

and  the  BCS  theory  of  superconductivity  clearly  indicate 

that  this  might  not  be  the  case.  As  a  matter  of  fact  some 

20 

of  the  physical  quantities  calculated  from  the  above  theories 
have  essential  singularities  for  vanishing  g  which  indicate 
that  one  would  not  be  able  to  obtain  these  results  via  the 
conventional  perturbation  theory.  One  would  like  to  know 
whether  the  non-analyticity  of  the  solution  in  term  of  g 
is  of  a  fundamental  nature  or  an  inherent  property  of  the 
Bogoliubov  treatment. 

As  an  example,  in  Chapter  2  we  construct  a  non-trivial 

many-body  model  Hamiltonian  for  which  the  S-matrix  can  be 

calculated  exactly.  We  then  investigate  the  analytic 

properties  of  S  as  a  function  of  the  coupling  constant  g 

and  the  number  of  particles  N  in  the  system.  In  Chapter  3 

we  construct  a  generalized  Bogoliubov  transformation  from 

21 

first  principles.  While  many  such  generalizations 


have 
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been  proposed  and  used,  none,  to  our  knowledge,  has  been 
derived  from  first  principles  and  none  is  as  general  as 
our  expansion.  It  is  instructive  to  show  that  such  an 
expansion  can  be  written  consistently  and  we  show  that 
the  usual  Bogoliubov  transformation  is  just  the  first 
term  of  this  infinite  expansion.  The  subsidiary  conditions 
that  result  from  this  expansion  are  discussed  in  Chapter  4. 

In  Chapters  5  and  6  we  solve  the  lowest  order  equation 
that  we  obtain  from  our  expansion.  In  Chapter  7  we  cons¬ 
truct  a  solvable  model  Hamiltonian  for  calculating  lifetimes 
of  quasi-particles  assuming  the  effective  interaction  between 
them  is  known.  In  Chapter  8  we  show  how  one  can  calculate 
this  effective  interaction  via  the  S-matrix.  In  Chapter  9 
we  define  a  new  boson  operator  expansion  which  leads  to 
what  we  call  the  G-matrix  equation.  With  this  approach  we 
can  calculate  improved  quasi-particles  energy  spectra  and 


their  lifetimes. 


. 


CHAPTER  2.  NON-ANALYTICITY  OF  THE  S-MATRIX 


IN  THE  MANY- BODY  SYSTEM 

For  a  given  interacting  system,  if  we  want  to  cal¬ 
culate  the  S-matrix  for  that  system,  using  perturbation 
theory,  then  we  must  assume  a  priori  the  analyticity  of 
the  S-matrix  in  the  coupling  constant  g,  for  vanishing  g. 

The  reason  for  this  is  due  to  the  following  three 
22 

theorems  : 

Theorem  I.  If  a  power  series  in  z  converges  for  some 
value  of  z,  say  z-^  ^  0 ,  then  it  converges  absolutely  for 
all  values  of  z  with  |z|  <  |  z^J  . 

Theorem  II ^  A  function  represented  by  a  power  series  is 
analytic  at  all  points  inside  of  its  circle  of  convergence 
Theorem  III.  A  function  analytic  at  a  point  may  be  repre¬ 
sented  by  a  power  series  expansion  about  this  point. 

Thus,  in  order  that  the  power  series 

oo 

s  =  l  Cn  gn  ,  (2.1) 

n=0 

that  we  obtain  from  perturbation  theory,  have  any  physical 
significance,  it  must  converge  for  at  least  one  value  of  g 
say  0.  Then  from  theorem  I  the  series  (2.1)  converges 

in  a  circle  centered  around  the  origin  with  radius  g^. 
Hence  from  theorem  II  we  conclude  that  S  is  analytic  at 
g  =  0.  If  we  calculate  the  S-matrix  exactly  and  find  that 
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it  is  not  analytic  around  g  =  0,  then  from  theorem  III 
we  see  that  we  cannot  represent  S  by  a  power  series  at 
this  point  and  hence  we  would  never  obtain  this  exact  S- 
matrix  from  perturbation  theory.  We  now  consider  just 
such  a  case. 


Consider  the  following  Hamiltonian 


N 

H  =  I  wi(aiai  +  +  g  I  <ai  +  al ) 2  (a-j  +  a^)2  Fii  (t) 

i=l  ii  j  j  j 


(2.2) 


where  (t)  is  as  yet  an  arbitrary  function  of  time.  Let 


1  ,  ^  .  ,  f1/  ,  9 

a_.=  —  (q-+  ip-)?  a.=  —  (q.-  ip.);  p.=  -i 

1  /2  3  3  3  /2  3  3  3 


3gj  • 


(2.3) 


Under  the  canonical  transformations  (2.3)  H  becomes 


H  =  H  +  Ht 
o  I 

where 


and 


(2.4) 


(2.5) 


(2.6) 


Thus  H  represents  N-harmonic  oscillators  interacting  via 
Hj.  We  want  now  to  calculate  the  matrix  elements  of  the 


10 


S-matrix.  For  simplicity  let  us  assume  that 


F±j  (t)  =  6  (t)  . 


(2.7) 


Recall  that  formally  (in  the  interaction  picture) 


S  =  P{exp[-i  /  H  (t')dt']} 


(2.8) 


where  P  is  the  chronological  time  operator.  In  our  case 


N 


2  2 


and 


S  =  exp  [-4ig  J  q.q.]  / 

•  •  JL  j 

il  J 


(2.9) 


S  =  <t  S  t  > 

mn  r  m 1  1  y  n 


(2.10) 


where  ijj  are  the  eigenfunctions  of  Hq.  Since  the  ground 
state  is  given  by; 


'f'o(ql-  •  •qN) 


N 

4  ,1,2^2 

=  77  exp  [-  2^1  +  q2  +  •• 


.qN)]  , 


(2.11) 


we  have  from  (2.10) 


N 

-  oo 

2  r  ,  r  ,  2  .  2  .  2 


Soo(q'N)  =  77  /  dqi  •  •  -  dqN  exp[-(q1  +  q?  +  .  ..qN) 


N 


L1  ^2 


N ' 


-  4ig(q^  +  q^  +  q^)2]  * 


(2.12) 


Thus 
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s 

oo 


(g,N=l) 


— i  /  dq  exp [-  (q2  +  4igq4)] 

/tt  0 

4  (iug)  2  exp [  jyjg] K\  ^32ig^  f°r  Im  g<0 


1  r  (~4i ) n  ( 4n )  !  n 

/~  n  n  i  g 

/tt  n=0 


(2.13) 

(2.14) 

(2.15) 


•  2  3 

where  K  is  the  Bessel  function  of  imaginary  argument 

Eq .  (2.15)  was  obtained  by  expanding  (2.13)  in  power  series 

of  g.  Similarly 


Soo(g'N=2) 


1  °°  2tt  2  4 

—  j  dr  f  d0r  exp[-r  -  4igr  ] 


TT 


(2.16) 


0  0 


=  (T6l?)"2eXp[I^I?]  [1  "  for  Im  g<  0 


v  (-4i ) n (2n) !  n 

=  l  - in - g 


n=0 


(2.17) 

(2.18) 


where  (2.16)  was  obtained  from  (2.12)  by  the  transformation 

q^=  r  sin  0  ,  q^  =  r  cos  0  ,  dq^dq^  rdrd0 .  $  is  the  proba- 

23 

bility  integral  .  As  can  be  seen  from  above  equations 
Soo  (g'N=1'2)  ps  smo°th  and  well  behaved  function  of  g,  for 
Im  g<0,  and  tends  to  unity  for  g ->  0 .  For  Im  g  >  0  it  diverges. 
The  radius  of  convergence  of  the  series  (2.15)  and  (2.18)  is 
zero  (ratio  test)  and  thus  SQO is  non-analytic  at  g  =  0  and 
it  cannot  be  obtained  from  perturbation  theory. 


Let  us  now  consider  Sqo  (g ,  N  >  3).  Define  the  gen- 

24 

alized  polar  coordinates  transformation 


12 


gl 

=  r 

sin 

i — 1 

CD 

sin 

0  2  ... 

sin 

CD 

** 

COS 

< P 

q2 

=  r 

sin 

i — 1 

CD 

sin 

Q  2 

sin 

9k 

sin 

<P 

q3 

_  r 

sin 

9i 

sin 

0  2  *  *  * 

cos 

9k 

(2.19) 

qk+1  =  r  sin  0^^  cos  02 


qk+2  =  r  cos  0 


where 


N  =  k+2  , 


k  =  1 , 2  ,  .  . 


(2.20) 


Then  from  (2.19)  we  can  show  that 


N 

l 

i=l 


=  r 


(2.21) 


Under  the  transformation  (2.19)  Eq.  (2.12)  becomes 


N 


2tt  tt 


7T 


Soo(g,N-3) 


=  TT 


f  dr  f  dcp  f  dd 


0 


0 


0 


,/dek  rk+l  X 


r  2  .  .  4 .  .  kQ  .  k-l_  .  _ 

exp[-r  -4i.gr  ]  sin  0-^sin  ©^...sin  0^ 
N  ,  i  k 

tt  ^  I(k)  J  exp[-x-4igx  ]  dx 
_  N  0  .  ,k  1 

TT  I  (k)  ( 8ig )  q  r(|+l)  X 


exp [32ig]  °_k _  ,  vr  8ig 
2 


(2.22) 

(2.23) 


)  for  Im  g  <  0  (2.24) 
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where  D  is  the  parabolic  cylinder  function  and 


-*  *  *  »  *  4 


13 


k  ^  k  —  1 

I(k)  =  J  sin  0  d0,  J  sin  0^d0 
0  1  1  0 


TT 


2  2  ' 


.  .  /  sin0,d0, 
0  K  k 


TT 


k 

2 


N 


-  1 


TT 


r<!  +  n 


r  (|) 


(2.25) 


(2.26) 


Hence 


Soo(g'N  *  3) 


N 


(8ig)  4  expt^]  D  N  >I8ig 


(/~i 


)  for  Im  g<  0 

(2.27) 


n 


v  (-4i) “  _  .  _  ,  N,  n 

l  '  -n."  '  F  (2n+  9  * 


r (|)  n=0  ni 


(2.28) 


Thus  S  (g,N  >  3)  is  smooth  and  well  behaved  function  of  g, 
for  Im  g<  0,  and  tends  to  unity  for  g ->  0 .  Here  again  the 
radius  of  convergence  of  the  series  (2.28)  is  zero.  Thus 
we  see  that  for  a  given  finite  N,  S(q,N)  is  a  non-analytic 
function  of  g  at  g  =  0 . 

From  the  definition  of  D  one  can  show  that 


D  N  [  (8ig)  ^2]  =  2  1} 

~2 


-1. 


(N-l)  [(16lg)  ]  (2.29) 


Isl^-oo 


!  ,5<N  X)  [16ig  (N-l)  ]  k  exp[~(N-l)-i(N-l)lni(N-X)-2(|7r|)!'2] 


N 


->  (16igN)  %  (|)  4 


(2.30) 


o ,  o 


Thus  as  N  00 ,  S 


->  0  (see  Appendix  A) 
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In  order  to  show  that  the  non-analyticity  of  S  is 
not  dependent  strictly  on  our  particular  choice  of  we 
have  considered  two  more  general  types  of  interaction: 


N 


2,  v 


1)  Consider  Hx  =  4g  (  V  q.)  6  (t)  with  N>2,  v  >  0. 

i=l 

Then  the  S-matrix  can  be  written  as 


°'°  r (|)  n=0  n! 


V  (~4i ) n  „  ,N  ,  ,  n 

l  — f  (T  +  nv)  g 


(2.31) 


The  radius  of  convergence  R  can  be  shown  to  be  (using  the 
ratio  test) 


,  .  T  (nv  +  ~  +  v) 

R  =  lim  -  - 


n->°° 


n 


TJ 


r  (Nv  +  j) 


n •  /  ,  N,-V,-l 

[lim  (nv  +  y)  ] 


n->oo 


(2.32) 


,  .  .  v  v-1 

lim  4v  n 

n->°° 


(2.33) 


Thus  for 


0  <  v  <  1;  R  =  °°,  Sqo  is  analytic  everywhere 


v  =  1  ;  R  =  i,  S  is  analytic  for  Igl  <  i 

'  4  o,o  J  11  4 


v  >  1  ;  R  =  0,  S  is  non-analytic  at  g  =  0  . 

0,0 


N 


2) 

Then 


Let  H-j.  =  4ig  l  q^  6  (t)  with  v  >  0 

i=l 


N 

oo  2  .  .  V  N 

S  =  ^  2  [  /  dqe^  “4l^  ]  . 

0/0 


(2.34) 


. 


> 
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The  radius  of  convergence  of  the  expression  inside  the 
bracket  (when  expressed  in  power  series  of  g)  can  be 
shown  to  be  (ratio  test) 

“  -  -  1 

R  1  =  lim  4 {j) 2  n2  .  (2.35) 

n->co 

Thus 


0  <  v  <  2  ,  R 

v  =  2  ,  R 

v  >  2  ,  R 


00 


1 

4 ' 


0, 


s 


0,0 


S 

0,0 


S 

0,0 


is  analytic  everywhere 
is  analytic  for  |g|  <  ^ 
is  non-analytic  at  g  =  0  . 


We  see  that  the  non-analyticity  of  the  S-matrix  in 
the  coupling  constant  is  a  property  of  many  systems (Appendix  B) . 
Although  a  rigorous  proof  does  not  exist,  there  are  indica¬ 
tions  that  a  similar  situation  prevails  in  most  non-trivial 
25  . 

theories  .  However  these  theories  might  still  be  analytic 
in  N  ^  or  in  the  product  gN  and  thus  in  term  of  this  constant 


we  might  still  be  able  to  employ  perturbation  methods. 
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CHAPTER  3.  QUASI-PARTICLE  EXPANSION  FOR 

A  BOSON-SYSTEM 


The  Hamiltonian  for  a  system  of  many  bosons  inter¬ 
acting  through  two-body  forces  is  given  by 


H  =  l  (£k'>j)akak+  1 


+  + 

„„  y  6  .  .  „  .  v(j-m)a.a'a  a 

2 ft  .  nL  ]  +  £,m+n  2  j  Z  m  n 

j fmn  t  ~  ~  ~  ^  ~  ~ 


(3.1) 


where  is  the  kinetic  energy  of  the  particle  with  momen- 

+ 

turn  k,  y  is  the  chemical  potential  and  a^  and  a^  are  the 
creation  and  annihilation  operators  satisfying  the  commu¬ 
tation  relation 


[ak,  at] 


=  6 


k,  j 


Also 


-ir.  (j-m)  ^  3 

v(j-m)  =  /  e  ~  ~  v(r)d  r 


(3.2) 


(3.3) 


If  one  assumes  v(r)  =  v(|r|)  then 

v  (q )  =  v(-q)  .  (3.4) 

t 

The  Heisenberg  equations  of  motion  for  a^  and  a^. 
can  be  derived  from  the  definition  of  the  time  derivative 
of  an  operator,  i.e. 


. 
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=  i [H,  ak]  ;  a£  =  i[H,  a^] 


t 


(3.5) 


From  (3.1)  and  (3.5)  we  find 


a  =  -i{  (e  -y)a  +  777-  Y  6  .  ,  n  ,  v(j-m)  x 
q  q  <3  jtan  ~  - 


<V  -V  <v  <v 


[6.  ata  a  +6.  ata  a  ]}. 
j  ,q  £  m  n  £  ,q  j  m  n 


“V 


With  the  help  of  (3.4)  we  find 


•  1  f 

a  =  -i{  (e  -y)a  +  r  I  6  ,  .  ,  v(q-m)aTa  a  }  . 

q  q  q  «  j^n  5+3'™+?  ~  ~  3  m  nJ 


<N/ 


<V  "V  <V 


(3.6) 


+ 

Suppose  we  define  A^  and  to  be  creation  and  anni- 

~  <v 

hilation  operators  such  that 


N 


(3.7) 


is  the  number  operator.  Furthermore,  assume  the  Hamiltonian 

■j* 

(3.1)  written  in  terms  of  A^  and  A^  commutes  with  N,  i.e. 

H  will  have  the  form 


t 


H  =  l  e,  a;  A,  + 


k  k  “k 


ft  ^  AkA£AkA£ek£  + 
k ,  £  ~  ~  ~  ~  ~~ 


(3.8) 


+  + 

If  one  finds  A^  and  A^  m  terms  of  a^  and  a^  and  determines 


the  coefficients  ^  . . .  etc.,  then  the  problem  is 


■ 


••  -  --  - 
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completely  solved.  Since  such  a  complete  diagonalization 
of  the  Hamiltonian,  except  for  some  solvable  models,  at 
present  seems  impractical,  one  can  try  to  find  approximate 
methods  of  diagonalization.  Among  these,  the  quasi-particle 
method  of  Bogoliubov  seems  to  be  the  simplest  one.  The 
method  considered  here  is  an  extension  of  a  technique  due 
to  Valatine"^ . 

+ 

We  define  b^  and  b^  as  the  creation  and  annihilation 
operators  for  quasi-particles  of  momentum  k  if  they  satisfy 
the  following  conditions: 


a)  -  Commutation  relation: 


b)  -  Time  derivative  of:  b^. 


“kbk 


-“kbk 


(3.9) 


(3.10) 


•  •  t 

c)  -  Certain  symmetry  properties  and  invariances  of  H,  a^ 

and  a^  should  also  be  satisfied  by  the  equations  expressing 
~  t 

these  operators  in  terms  of  b^  and  b^  (see  below) .  In  this 
way  the  quasi-particle  operators  are  characterized  by  one 
parameter,  which  is  associated  with  a  quasi-particle 

in  the  state  k.  The  parameter  co^  will  play  the  same  role  as 
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in  Eq.  (3.8).  For  exactly  solvable  problems  however,  the 
equation  of  motion  for  is  not  the  same  as  that  of  b^ 
i.e.,  instead  of  (3.10)  we  have 


[  e  + 
q 


-  y 
« 1 


% 


'ql 


(AIA£}  + - ]Aq  * 

'V  'V 


(3.11) 


Let  us  consider  a  system  of  interacting  bosons. 

The  bulk  properties  of  such  a  system  (like  pressure)  depend 
on  the  energy  level  spectrum  of  the  particles.  Therefore 
we  can  conceive  of  a  model  of  this  system,  where  different 
states  (quasi-particles)  are  non-interacting  and  are  charac¬ 
terized  by  their  energies  co^, .  This  model,  which  is  useful 
in  calculating  the  gross  properties  of  the  original  system, 
has  a  relatively  simple  mathematical  structure.  However, 
it  will  not  provide  a  complete  description  of  the  actual 
system,  because  there  will  be  interaction  between  the  quasi¬ 
particles  caused  by  the  residual  part  of  the  Hamiltonian. 
Thus,  we  want  to  find  the  operators  b^  such  that  will  be 

■V  *■>* 

very  close  to  the  actual  energies  of  the  different  states, 
and  then  consider  the  weak  interaction  between  the  quasi¬ 
particles.  To  this  end  we  assume  that  as  t  -*  ±<»  we  have  a 

system  of  non-interacting  quasi-particles,  and,  in  this 

+ 

limit,  the  operators  a^  and  a^.  can  be  expressed  as  linear 

+ 

combinations  of  b,  and  b,  i.e. 
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lim 

t->±CO 


ak  — ►  h(k)bk 


+  g(k)b+, 
~  -k 


t 

lim  a, 

.  k 


h(k)b£  + 


g  (k)b 


-k 


(3.12) 


where  h(k)  and  g(k)  are  assumed  to  be  real  and  symmetric 
coefficients  to  be  determined.  The  transformations  (3.12) 
for  k  ^  0  take  into  account  the  fact  that  the  propagation 
of  a  single  particle  by  itself,  in  the  medium,  is  no  longer  mean¬ 
ingful,  due  to  the  interaction,  and  there  is  a  back  flow  of  other 
particles  around  it  as  it  moves  through  them.  If  a  particle  has 

no  motion  (k  =  0) ,  there  is  no  reason  why  a  should  have  the 
~  o 

same  relation  to  b  as  for  other  values  of  k.  However,  it 

o  ~ 

is  convenient  to  assume  that  Eqs .  (3.12)  are  valid  for  all  k. 

We  shall  see  later  on  that  this  will  mean  that  we  have  to 
impose  a  subsidiary  condition  on  g(0). 


To  find  the  complete  expansion  of  a^  for  finite 

values  of  t,  we  must  first  study  the  symmetries  and  the 

invariances  of  such  an  expansion.  Now,  if  by  a  unitary 

+ 

transformation  of  a^  and  a^ ,  which  does  not  depend  expli¬ 
citly  on  time,  the  Hamiltonian  (3.1)  or  the  equation  of 
motion  (3.6)  remain  invariant,  then  the  invariance  should 

be  preserved  as  t  ->  ±°°  .  This  implies  that  the  Hamiltonian, 

+ 

expressed  in  terms  of  b^  and  b^.  must  have  symmetries  cor¬ 
responding  to  the  unitary  transformation  under  consideration. 


' 

- 


21 


Let  elS  be  the  generator  of  such  a  transformation,  then  H 
must  commute  with  e1^,  i.e., 


[H,  exp  (is  (a,a^,k)  )  ]  =  0  (3.13) 

i 

Here  S(a,a,k)  denotes  that  S  is  a  functional  of  all  of  the 

+ 

operators  a^  and  a^  and  in  addition,  it  may  also  depend  on 
the  momenta  of  the  particles  in  the  system. 


The  following  two  symmetry  properties  of  the  Hamil¬ 
tonian  are  useful  in  deciding  the  form  of  the  expansion. 

1)  -  Galilean  transformation  -  Let  S  be  defined  in  such  a 
way  that 


e 


iS 


-iS 

a  e 

q 


imv .  q 

a  e 

q 


e 


iS 


- imv . q 


(3.14) 


where  v  is  a  constant  (velocity)  vector  and  m  is  the  mass 
of  a  particle.  The  Hamiltonian  (3.1)  and  the  equation  of 
motion  (3.6)  remain  invariant  under  this  transformation. 
An  explicit  form  of  the  operator  S  can  be  obtained  by  the 
method  outlined  below.  Here  we  give  the  result: 


+  + 

=  -m [  l  (k.v)a  a  -  £  (k.v)a_  a_  ] 

k>0  ~  ~  5  ~  k>  0  ~  ~  £  £ 


S 


(3.15) 


' 


. 
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2)  -  Invariance  under  a  unitary  transformation  elU  such 
that 


-iU 

e 


(3.16) 


(3.17) 


Again  H  and  the  equation  of  motion  (3.6)  remain  invariant 
under  this  transformation.  A  method  for  constructing  the 
operator  U  is  as  follows:  We  write  equation  (3.16)  as 


i  (2n+l )  tt 

tt  d 

q 


n=0 , ±1 , ±2 ,  . . .  . 


Then  equation  (3.17)  has  the  form 


(3.18) 


e 


iU 


e 


-iU 


-i  (2n+l)  tt  t 
0  0 

q 


(3.19) 


i.e.,  once  we  fix  the  phase  in  equation  (3.16)  the  phase 
in  equation  (3.17)  is  fixed.  Thus  equations  (3.18)  and 
(3.19)  are  consistent.  Now  rewrite  (3.18)  as 


iU 

e 


e 


109  a9 


log  aq  i  (2n+l)  tt  iU 
e  -?  e  e 


(3.20) 


If  we  assume  that  the  commutator  of  U  and  log  a  is  a  c- 

q 


number  then  one  has 
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a 

q 


log  aQ  iU  i[U'lo9 
e  ~  e  e 


(3.21) 


Comparing  (3.21)  and  (3.20)  we  find  that 


[U ,  log  a  ]  =  (2n+l)  tt  .  (3.22) 

q 

In  a  similar  fashion  one  obtains  from  Eq.  (3.19) 

[U ,  log  a+]  =  -  ( 2n+l )  tt  .  (3.23) 

q 

( 2  n  "i"  1 ) 

From  Eq.  (3.22)  one  sees  that  the  two  operators  - - - -  U 

TT 

and  log  a^  are  quantum  mechanical  conjugates  of  each  other. 
There  are  an  infinite  number  of  operators  conjugate  to 
log  a^  (corresponding  to  different  values  of  n) .  Thus 

U  =  {U  /[U  ,  log  a  ]  =  (2n+l)iT}  . 
n'  n  ^  q 

The  general  form  of  Un  is 

U  =  -  (2n+l )  it  [  l  a^a  -  l  a_  ]  +  F  (a.)  (3.24) 

n  k>0  ?  5  k>0  ~  5  n  5 

where  F  (a^)  is  an  arbitrary  function  of  its  argument. 

+ 

Since  from  Eqs.  (3.16)  and  (3.17)  one  can  show  that  U  =  U 

we  find  that  F  (a,  )  must  be  a  c-number  (we  choose  it  to 

n  k. 


be  zero).  Hence, 


' 
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U  =  -  (2n+l )  tt  [  l  a+a  -  J  a+  a  J  (3.25) 

n  k>0  ?  5  k>0  ~  5 

satisfies  (3.18)  and  (3.19)  as  can  easily  be  checked  by 
using  the  identity 


eXA  Be  XA  =  B  +  A[A,B]  +  ~  [A,[A,B]]  +  ...  (3.26) 


Since  U  commutes  with  the  Hamiltonian,  then  as  t  tends  to 
±°o,  the  invariance  under  the  change  of  sign  of  the  crea¬ 
tion  and  annihilation  operators  should  hold  even  when  H 

J. 

is  expressed  in  terms  of  b^  and  b^ .  To  prove  this  we  want 
to  show  that 


lim 

t-*±o° 


e 


-iU 


lim 

t->±00 


e 


-iU 


(3.27) 


We  note  that  Eqs.  (3.27)  imply  that  U  should 
functional  form  (up  to  a  c-number)  under  the 
(3.12).  Now  in  Eq.  (3.25)  the  term 


y  a,a,  -  J  a  ,a  , 
i  L _  k  k  .  £  -k  -k 
k>0  ~  ~  k>0 


have  the  same 
transformations 


transforms  as  required  (as  we  shall  show  below).  However, 

the  term  for  k  =  0  i.e.  a^a  does  not.  Thus  we  cannot 

o  o 

satisfy  Eqs.  (3.27)  for  q  =  0.  However,  Eqs.  (3.27)  are 
valid  for  q  ^  0.  To  see  this,  substitute  Eqs.  (3.12)  into 
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into  Eq .  (3.25)  and  obtain 


U  (t  -*  ±°° )  — >-(2n+l)7i{  V  [h2  (k)-g2  (k)  ]  (±«>)b,  (±<*>) 

k>  0  ~  £  5 

-  I  [h2(k)-g2(k)]b  +  ,  (±~)b  ,  (±«)  +  F(b  ,b+)}  (3.28) 

k>  0  15  "C  °  ° 

where 


F(b  ,b+)  =  h2(0)b+b  +  g2(0)b  b++  h(0)g(0)  (b  b  +  b+b+) 
o  o  oo  3  o  o  ^  o  o  oo 


Eq.  (3.28)  can  be  simplified  if  we  note  that 


(3.29) 


h2  (k)  -  g2  (k)  =  1  ,  (3.30) 

which  is  the  condition  that  the  transformation  (3.12)  be 
a  canonical  transformation,  i.e.,  the  condition  that  both 
Eqs.  (3.2)  and  (3.9)  are  satisfied  when  t  ->  ±°°. 

t 

From  the  equation  of  motion  for  b^  and  b^  (Eq.  (3.10)) 

we  have 

bk(t)  =  bk(±T)  exp[-imk(t  +  T)  ]  , 

T-arbitrary  (3.31) 

bk(t)  =  bk(±T)  exp[iwk(t  +  T)  ]  , 
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from  which  it  follows  that 

bk(t)bk(t)  =  bk  (±T )  bk  (±T )  •  (3.32) 

Hence  we  can  write  Eq.  (3.28)  as 

U(t)  =  -  (2n+l )  tt  l  bUt)b,(t)  -  h*  (t)b.  (t) 

k>0  E  5  5  5 

+  F(bo(t)  ,b£(t)  )  (3.33) 

2  2 

provided  the  relation  h  (k)  -  g  (k)  =  1  can  be  satisfied 

for  all  values  of  t.  Note  that  we  used  in  (3.33)  the  fact, 

which  we  shall  discuss  in  a  later  chapter,  that  go  =0. 

o 

The  validity  of  (3.30)  for  all  times,  will  be  imposed  on 
the  expansion  of  a^.(t)  as  a  subsidiary  condition.  Thus 
we  can  write 


iU(t) 


bq(t) 


-iU(t) 


b  ( t) 

q 


q  t  0 


(3.34) 


Similarly  we  also  have 


is  -is  imY‘S 

elb  b  (t)  e  lb  =  bg(t)  e  V  q  (3.35) 

(Note  that  Eq.  (3.35)  is  true  for  all  q  because  of  the 
term  k.v  in  Eq.  (3.15)).  Since  H  commutes  with  U  and  S, 
the  transformations  (3.34)  and  (3.35)  should  leave  H 


»  •  •>.  *  -  .  -It 
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invariant.  This  means  that  in  the  expansion  of  H  in 

•j* 

terms  of  ' s  and  bk '  s  every  term  must  contain  an  even 

number  of  creation  and/or  annihilation  operators  (a 

consequence  of  Eq.  (3.34)).  The  correct  dependence  upon 

k  will  be  deduced  from  Eq.  (3.35).  The  same  argument  can 

+ 

be  applied  to  the  expansion  of  a^  and  ak ,  with  the  result 
that  each  term  in  the  expansion  has  an  odd  number  of  b's. 
The  most  general  expansion  which  preserves  the  above 
mentioned  invariances  to  each  order  is 


ak  =  h(k)bk  +  g(k)b^k  +  e  a  1  ^ 1  ( — ^—3-)^  /  dJxd°yd 


-a  t 


ft  %  2 


,3  , 3  -,3 


(2tt  ) 


x  {f  (x^/Zjb^.b  b  +  f  9  (x,y,  z)b+_b+  b 
1  ~  ~  ~  -x.  v  z  z  ~  ~  ~  ~x  -y  z 


+  f3(5'y'5)bxbybz  +  f4(?'^!)b-xb-yb-z} 


x  6 (k-x-y-z)  +  terms  having  five  operators  +  ...  (3.36) 


where  ft  is  the  volume  of  quantization  and  a  is  a  small  (a<<^k) 

positive  number.  The  functions  h(k),  g(k),  f^,  f2,  f3 

and  f .  are  all  assumed  to  be  real  functions  of  their  argu- 

ments.  We  shall  assume  that  this  expansion  is  valid  for 

all  k  including  k  =  0  (see  paragraph  after  Eq .  (3.12)). 

-j*  , 

From  (3.36)  we  find  the  expansion  for  ak  ,  which  will 
depend  on  the  same  functions  as  ak  does.  Equation  (3.36) 
together  with  its  complex  conjugate,  and  the  definition 


of  the  time  derivative  of  b^.  (Eq.  (3.10))  will  give  the 

~  i 

complete  expansion  of  the  boson  operators  a^.  and  a^.  in 
terms  of  the  quasi-particle  operators  b^.  and  b^. .  The 
complete  dynamics  of  the  system  can  be  specified  by  the 
equation  of  motion  (3.6),  the  canonical  commutation 
relations  (3.2)  and  (3.9)  and  by  the  expansion  (3.36). 
The  physical  significance  of  the  functions  h,  g  and  the 


f's  will  be  considered  later. 


\ 
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CHAPTER  4.  COEFFICIENTS  OF  THE  EXPANSION 

If  we  truncate  the  expansion  of  a^  (Eq.  3.36)  and 
keep  all  of  the  terms  containing  products  of  m  operators 
or  less,  then  the  number  of  the  coefficients  of  the  ex¬ 
pansion  is  ^(m+l) (m+3)  (m=l , 3 , 5 . . . ) .  The  non-linear 
integral  equations  satisfied  by  these  coefficients  are 
determined  by  substituting  Eq.  (3.36)  and  its  conjugate 
in  the  equation  of  motion  (3.6)  and  arranging  the  products 
of  the  operators  in  the  normal  order.  By  equating  the 
coefficients  of  the  same  set  of  operators  on  the  two  sides, 
after  using  (3.10),  we  find  -^-(m+l)  (m+3)  coupled  integral 
equations  for  the  unknown  functions.  Since  the  transfor¬ 
mation  of  a^  to  b^  is  a  canonical  transformation,  the 
commutation  relation  (3.2)  should  be  preserved  to  the  same 
order  of  the  expansion.  This  means  that  on  •^•(m+1)  (m+3) 
functions  we  have  a  number  of  subsidiary  conditions.  For 
example  if  m=l,  we  have  two  functions  h(k)  and  g(k)  and 
one  subsidiary  condition  (3.30).  In  general  by  keeping  m 

operators  in  the  expansion,  we  find  from  Eqs.  (3.2)  and 

12  .  . 

(3.9)  that  there  are  ^-(m+l)  of  these  conditions. 

Thus  by  substituting  (3.36)  in  (3.6),  eliminating 
the  time  derivatives  of  bk  using  (3.10)  and  then  equating 
the  coefficients  of  the  normal  products  of  b^ ' s  we  find 
for  m=3  (finite  time) 


. 

;■ 

. 

’  -  '  -V.  «  -#*  «  —  -*  ■*  «  *  -  —  - 

. 

°  ,t<'  I 
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[a)k"£k+p]h(k)  =  I  [v  (0  )  +v  (k— j)  ]g2  (j)  }h(k) 

~  ~  j 


+  I  v  (k- j  )  h  ( j  )  g  ( j  )  }g(k) 

j 


+  ^  ^k+ j ,m+n 

jmn  ~  ^ '  ~  ~ 


v (k-m) { 2f „ ( j  ,-m,k)  g  ( j  )  h  (m) 


+  2f i (3 r-n^k) g ( j ) g (n) +6f 4 (k ^n^m) g (m) g (n) }  + 


(4.1) 


^Wk  +  ek_y  g  ^  =  I  [v  (0) +v  (k- j  )  ]  g2  ( j  )  }g2  (k) 

~  ~  ~  j 


“  (§  I  v  (k- j  )  h  ( j  )  g  ( j  )  }h(k) 

j 


h  J  6q+j/m+n  V(k"?){2f2(k'2'"^)g(2)g(?} 
jmn  ~  ^ ' ~  ~ 


+  2f2 (n,m,-k)g (m)g(n)+6f4 (j ,-m,k)g(j)h(m) }  + 


[o3x+Wy-toz-ex+y+z+y]  f±  (x,y,z)  =  ^  {  v  (k-y)  h  (x)  h  (y)  h  (z ) 


(4.2) 


+  v (k-z ) g (x)g (y)g  (z)  +  v (k-x) g (x) h  (y) g  (z ) }  + 


(4.3) 


[03x+a)y-a3z+ex+y+z-y]  f2  (x,y,z)  =--  {  v  (k-y )  h  (x)  g  (y )  h  ( z ) 


+  v(k-z)h(x)g(y)h(z)  +  v  (k-x)  g  (x)  g  (y )  g  ( z )  + 


(4.4) 


' 
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^x+(x)y+toz-ex+y+z+VJ^  f2  Q'X'  ^  =  §  v  (k-y)  g  h  ^  h  ^  +  *** 

<v  -v  <v  -v  <v 

(4.5) 

[wx+a3y+cOz+ex+y+z“Vi]  f4  (x,y,z)  =  -i  v (k-y)h (x) g  (y) g  (z)  +  ... 

A#  <V  'V 

(4.6) 

The  number  of  terms  in  these  relations  (especially  for  f ^ , 
f^  and  f^)  is  large  and  we  have  given  only  the  leading 
terms  in  each  case. 


There  are  four  subsidiary  conditions  on  h(k),  g(k),  f^. 


■F  F  *F  • 

r  2  r  3  '  r  4  * 


1.  h2  (k)  -  g2  (k)  =  1 


(3.30) 


2. 


f3(x,y,z)  -  f2(x,y,z)  =  0 


(4.7) 


3.  f1  (r,s,k)h  (k)  -  f  2  (k  ,r ,  s)  g  (k)  +  f ± (-s ,-r ,q) h (q) 


-  f  (q,-s,-r)g(q)  +  — — ^  fd3yd3z  6  (k+s-y-z) 
^  ~  ~  ~  ~  (2tt) 


[f1(-s,y,z)f1 (r,y,z)  -  9f 4 (y , z ,-s) f 4 (y , z ,r) ] 


/ d3yd3z  6  (k+r-y-z)  [f  (y , z ,-r) 


(2tt) 


x  f  (y ,  z ,  s)-9f3  (y,z,-r)f3 (y,zrs) ] 


=  0  ;  q-k  =  r+s 


(4.8) 


' 
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4. 


+  3f3  '“r)h(g) 


3  3 

/ d  yd  z  6  (k+s-y-z) f 3  (y, z ,-s) f 1 (r,y, z) 


(^tt  ; 


0 


q-k  =  r+s 


(4.9) 


The  eigenvalues  can  be  found  from  Eqs.  (4.1)- (4.6)  by 
first  solving  f^,  f ^ ,  f3  and  f^  in  terms  of  h(k)  and  g(k). 


equations,  and  therefore  their  solutions  will  be  func¬ 
tionals  of  h(k)  and  g(k).  By  substituting  f^,  f^,  f^  and 
f^  in  Eqs.  (4.1)  and  (4.2)  we  find  a  coupled  set  of  non¬ 
linear  equations  which  determines  h(k)  and  g(k)  and  the 
eigenvalues  oo^..  Except  for  the  lowest  order  term,  i.e. 
for  m  =  1,  it  is  difficult  to  consider  the  question  of  the 
compatibility  of  the  solutions  with  the  subsidiary  condi¬ 
tions  obtained  from  the  requirement  that  the  a's  and  the 
b's  both  satisfy  canonical  commutation  relations.  In 
practice,  it  may  be  more  convenient  to  find  approximate 
solutions  of  Eqs.  (4.3)- (4.6)  which  are  compatible  with 
these  subsidiary  conditions. 
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CHAPTER  5.  QUASI-PARTICLE  ENERGIES  IN  THE 

LOWEST  ORDER  APPROXIMATION 


In  the  lowest  order  of  approximation  (m  =  1)  we 
have  the  following  equations  for  h(k)  and  g(k) 


[^k“  £k+  y  -  A  (k)  ] h (k)  =  A (k) g  (k) 


(5.1) 


[o)k+  ek"  y  +  A  (k)  ]  g  (k)  =— A  (k)  h  (k)  , 

where 

A(k)  =  i  J  [V  (0)  +  V(k-S>)]g2U) 

“  £ 

and 

A (k)  =  b  l  V (k-£)h (£)g (£) 

“  0. 


(5.2) 

(5.3) 

(5.4) 


Equations  (5.1)  and  (5.2)  are  obtained  from  Eqs.  (4.1) 
and  (4.2)  by  neglecting  terms  containing  f^,  f^r  f-^  an<3 
f^.  The  linearized  form  of  the  equations  (5.1)  and  (5.2) 
yield  the  energy  spectrum  as 


“k  =  {[£k  _  y  +  A(t)]2  ■  t2(V}h 


(5.5) 


Note  that  equation  (5.5)  does  not  give  ook  explicitly  since 
A(k)  and  A(k)  depend  on  Multiplying  (5.1)  by  g  (k) 

and  (5.2)  by  h(k)  and  adding  the  resulting  equations  we 
find,  using  Eq.  (3.30): 


. 

.  ' 
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A  (k) 

g  (k)  h  (k)  =  -  .  (5.6) 

k 

From  (5.4)  and  (5.6)  we  obtain  an  integral  equation  for 
A  (k) 


A  (k) 


_1_ 

2ft 


J  v(k-£) 


A(£) 


(5.7) 


The  solution  of  this  equation  is  subject  to  the  subsidiary 
condition  (3.30).  In  addition,  the  number  of  particles 
in  the  system  is  given  by  the  expectation  value  of  the 
number  operator. 


n 


a+a.  I  0> 
1  t  K  K  1 

k  ~  ~ 


o  i  y 


n  +  <  0 
o 


y 

k^O 


a.^a,  I  0>  . 
k  k  1 


(5.8) 


where  | 0>  denotes  the  ground  state  of  the  entire  system, 
and  is  the  vacuum  state  for  the  quasi-particle  operator. 

•4*  , 

Using  the  lowest  order  expansion  of  a^  and  a^  we  find 
that 


n 


n 


+  I 

k^O 


g2  (k) 


(5.9) 


Let  us  define 


n=g^(0).  (5.10) 

o  ^ 

Then 

n  =  l  g2  (k)  .  (5.11) 

k 


■ 
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Thus  g  (k)  is  the  average  number  of  bosons  in  the  state  of 
momentum  k. 

2  g 

Since  from  general  considerations  one  has  a)  =  0, 

we  must  treat  the  k  =  0  state  differently.  Thus  we  impose 

on  our  equations  the  condition  (5.11)  and  the  requirement 

goq=  0.  From  a>o=  0  we  determine  the  chemical  potential  y 

2 

and  from  (5.11)  we  determine  g  (0)  (=nQ) .  Note  that  since 

we  impose  the  condition  coq=  0,  we  must  be  careful  in  deal¬ 
ing  with  the  k  =  0  state  (for  example,  Eq.  (5.7)  is  valid 
only  for  k  ^  0).  Similarly,  as  we  shall  see  later,  Eq. 
(5.7)  must  be  modified  if  it  is  to  be  true  for  all  k. 

If  ft  is  the  volume  of  the  system  then 

F  (k)  =  S-g2(k)  (5.12) 

represents  the  normalized  distribution  function  for  bosons. 
Thus , 


- /  F(k)d3k  =  1 

(2rr) 

From  Eq.  (3.30)  we  obtain 


(5.13) 


=  [1  +  2.  F  (k)  ]  ^  . 


h  (k) 


(5.14) 


. 
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We  now  solve  the  coupled  eigenvalue  equations  (5.1)  and 
(5.2)  by  approximating  A(k)  and  A(k)  using  a  known  dis¬ 
tribution  function  F(k).  This  approximation,  in  effect, 
linearizes  the  equations  for  h(k)  and  g(k).  Thus  we  write 


A(k)  = 


A(k)  = 


1  r 

Replacing  —  l 

u  k 


£  [V(°)  +  ~  l  V(k-^)F(£)]  , 


~  l  V(k-£)  [F  01)  (§■  +  FU))]1"2 

by  — 1 — _  Jd3k,  we  find 
(2tt) 


(5.15) 

(5.16) 


A  (k)  =  ^  [V(0)  +  V (k) ]  ,  (5.17) 

A (k)  =  ^  V (k)  ,  (5.18) 


where 

V  (k)  =  — T  I  V  (k-Z)F  (Z)d3Z  .  (5.19) 

(2tt)J 

In  arriving  at  (5.18)  we  have  assumed  that  —  <<  F(£)  for 
all  values  of  Z.  Using  the  above  equations  and  Eq.  (5.5), 
we  obtain: 


{  le y  +  2.  V  (0)  +  2.  V(k)]2  - 


0) 


k 


n 


V2(k)}% 


(5.20) 


for 


wo= 


all 


k.  We  now  determine 
Since  sq=  0,  we  find 


the 


chemical  potential  y  from 


0. 


(5.21) 


y  =  +  £  v(0) 

and 

2 

Wk=  { [£k  +  ^  V(k)]2  "  ^2  V2  {~)}h 
If  we  choose  F (k)  to  be 


F  (k)  =  ( 2 tt )  3  6  (k)  , 

then  from  (5.19) 

V  (k)  =  V  (k)  , 

'V  <v 

and  (5.22)  now  becomes 

2 

wk  =  {  [£k  +  ^  V(t}  ]  2  "  ^2  y2  (t}  }1"2 
~  ^"2 


which  is  just  the  result  obtained  by  Bogoliubov. 


8 


(5.22) 


(5.23) 


(5.24) 


(5.25) 


CHAPTER  6.  SEPARABLE  INTERACTIONS 


We  shall  now  proceed  to  study  the  exact  solution 
of  Eqs.  (5.1)  and  (5.2).  However,  we  first  must  modify 
the  integral  equation  (5.7).  As  we  have  seen,  the 
requirement  that  go  =  0,  which  we  imposed  on  our  equations, 
implies  that  Eq.  (5.6)  is  valid  only  for  k  ^  0 .  However 
Eq.  (5.4)  is  valid  for  all  k.  Thus  we  rewrite  (5.7)  as 

<v 


A  (k)  =  i  V  (k)  h  (0 )  g  (0 ) 


2  ft 


V (k-£ ) 


A  U) 


(6.1) 


where  the  prime  on  the  summation  indicates  that  £  =  0  is 
omitted . 


Now  let  us  consider  a  non-separable  potential. 

Let  us  assume  that  it  is  square  integrable  (almost  all 
potentials  satisfy  this  summation  )  i.e. 

00  OO 

f  /  |  <  k  |  v  |  q>  |  d  k  d  q  <  00  .  (6.2) 

0  0  ~ 

Define  a  set  of  functions  g  (k)  by  the  integral  equation: 

°°  3 

gn (k)  =  An  /  1 vlg>gn(9)  d  q  *  (6,3) 

27 


Then  according  to  Mercer's  theorem 


■ 


■ 


Since  this  series  is  convergent  we  can,  therefore,  appro¬ 
ximate  any  non-separable  potential,  very  well,  by  a  finite 
number  of  separable  terms.  Hence  let 
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n 

V (k-£  )  =  Y  A  a  (k )  a  (£ ) 

~~  L  nn~  n~ 

n=l 

=  A^a-^  (k)ot^  (£ )  +  <f>^(k,£)  , 

where  n'  is  a  given  number.  Substituting  Eq.  (6 
Eq.  (6.1)  we  obtain 

1  A  (£) 

A(k)  =  f(k)  -  ciai(k)  -  2si  r  — 

X/  X/ 

where 

f  (k)  =  —  [A-^a-^  (k)a^  (0)  +  4)-^  (k,0)  ]h  (0)g  (0) 

and 

An  a  (£ )  A  (£  ) 

-  V  1  1  ~ 

C1  ~  l  oo 

£  ~ 

Define 

Aj^  (k)  =  A  (k)  +  c1a-L  (k) 

Then  from  (6.8)  and  (6.9)  we  have 
An  An  ,  a?(k)  - 

ci=  d  11  +  2?r  I 

From  Eqs.  (6.6),  (6.9) 

equation  for  A-^(k); 


V 

£ 


a1  (£)  A-l  (£) 


0) 


£ 


(6.5) 

.5)  into 

£ )  ,  (6.6) 

,  (6.7) 

(6.8) 

(6.9) 

(6.10) 


and  (6.10)  we  find  an  integral 
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1  A ,  (£) 

A  (k)  =  f  (k)  -1  ' 

1  ~  ~  p  ^ 

Ax  ~  ,  a  1(A')<f)  1  (kr A 1 ) 


(k ,  £)  - 


2ft  al  £ 


co 


ahr) 

i  +  4  r  1  - 


A. 

2ft 


£'  W£ 


(6.11) 


We  now  write 


<j>^(k,£)  =  A2a2(k)a2(£)  + 


(6.12) 


and  repeat  the  whole  procedure  till  we  run  through  all 
the  terms  in  (6.5).  Thus  if  n'=  2  then  (f>2(k,£)  =  0  and 
equation  (6.11)  becomes 


A^(k)  =  f  (k)  +  c2a2  (k)  , 


(6.13) 


where 


= 


,  A,  (£) 

-  I  ■  1,,~-  D  ( £ )  , 


£ 


co 


£ 


(6.14) 


and 


D(£)  = 


h 

ot  ^  ( £ ) 

a  n  (£ )  a0  (£  '  )  -i 
^  •  jl  ~  z  ~ 

A2 

ot  2  ( £ ) 

2ft 

£■  “i- 

2ft 

1 

A1 

.  1 

a?(£') 

j* '  1  ~ 

2ft 

V  aJ£'  J 

(6.15) 


Substituting  (6.13)  into  (6.14)  we  have 

f  (£) 


C2 


-  r  u. 
£  £ 


D  (£) 


i  +  r 

i 


a?  (£ ) 


co 


£ 


(6.16) 


D  (£) 
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With  the  help  of  Eqs.  (6.9)  and  (6.13)  we  finally  have 


A (k)  =  f  (k)  -  c.a,  (k)  +  c0a0(k)  . 

~  JL  _L  ~  Z  Z  ~ 


(6.17) 


Let  us  now  consider  the  simplest  form  of  a  separa¬ 
ble  interaction; 


v (k-£)  =  Xa  (k)a  U) 


(6.18) 


Then  the  solution  of  Eq.  (6.1)  can  easily  be  found  to  be- 


I  “<*>  ro 


A  (k)  = 


\  r*  * 

1  +  ^  z 


a2  (£ ) 


d2  a  (k)  , 


Z  ^Z 


where 


(6.19) 


T  =  h  (0)  g  (0)  a(0)  , 

o 


(6.20) 


and 


do  = 


—  r 

n  o 


i  +  _L  y 

1  90  i 


a2  (£) 


2fi  Z  W Z 


(6.21) 


Also  we  note  that  Eq.  (5.3)  can  be  written  as 


A(k)  =  a2(0)  +  d^  a  (k)  , 


(6.22) 


where 


d.  =  £  I  aU)g2U) 
1  u  all  Z 


(6.23) 


' 
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and 

n  =  I  g2  U)  •  (6.24) 

all  l 

From  Eq.  (5.5)  we  have 

“k  =  {  [£k'  p  +  JT  “2  (°)  +  ^ct  (k)  ]  2 
Imposing  the  requirement  c oq=  0,  we  get 

y  =  ~  a 2 ( 0 )  -  (d2-  dx)  a ( 0 )  .  (6.26) 

Hence 

2  2  2  t 

Wk=  ^£k+  (d2-dl^ a  ^  +  dla  ^  ~  d2a  ^ 2  (6.27) 


.,2  2  n  .  -i 

-  d2a  (k) } 

(6.25) 


In  order  to  determine  the  distribution  function  F(k)  we 
eliminate  h(k)  between  Eqs.  (3.30)  and  (5.6)  to  obtain 


4  2  a2»> 

g  (k)  +  g  (k)  -  - —  =  0  ,  k  ^  0  .  (6.28) 

4co, 

k 

This  quadratic  equation  has  the  acceptable  solution 


g2(k) 


,  A2  (k)  % 

j  {  [X  +  - 2~]  ”  1}  k  ?  0 

“k 

e  -  y  +  A  (k) 

i  — v — -  - 1}  *  * 0 


(6.29) 
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We  then  determine  g2(0)  (and  hence  h2(0)  since  h2(0)  = 

1  +  g2  (0 ) )  from 

g2 (0)  =  n  -  l '  g2  (k) 
k 

=  n  -  l  ±{[1  +  — j — ]  -  1}  .  (6.30) 

k  uk 

2 

Let  us  discuss  the  dependence  of  g  (k)  on  k.  For  a  short 
range  potential 

0  - 

which  implies  that  tends  to 

k2  ~ 
ek  ~  2m 

for  large  values  of  k  (see  Eq.  (6.25)).  In  this  limit 
(see  Eq .  (6.19)) 


A2  (k) 

— f—  «  1  ,  (6.31) 

wk 


and  by  expanding  Eq .  (6.29)  we  find  the  asymptotic  form 

of  g(k) 


g2  (k) 


1  <k) 

k->°°  4  2 

0). 


m  d2  a2 (k) 


(6.32) 


For  a  potential  V(r)  which  depends  only  on  the  distance 
of  the  particles ,  the  Fourier  transform  V(k)  is  a  function 


. 
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of  k  only.  It  can  be  obtained  either  as  a  power  series 
2 

in  k  ,  or  as  a  limit  of  such  V(k) 's.  Thus  we  shall  assume 

2 

that  our  a (k)  is  a  function  of  k  only.  Hence  as  k  goes 
to  zero  by  using  a  Taylor  expansion  of  the  expression 
inside  the  square  root  of  Eq.  (6.27)  we  obtain 


co 


k  k+ 


+  (d  -  d  )  (^2L.)  a(0)]%  (2d  )55.  (6.33) 

zm  l  z  dkz  kz=Q 


Rewriting  Eq.  (6.29)  as 


2  / 1  \  1 
9  <*>  =  25- 


{  [-2  +  A2  (k)  J  ^  -  wkJ 


(6.34) 


we  have  thai 


g2  (k) 


|d0a(k) 


k->0 


to. 


(6.25) 


2 

provided  a(0)  f  0.  Note  that  in  order  that  g  (k)  be 
finite  for  k ->  0  we  must  treat  the  k=0  state  differently. 
We  also  note  that  for  small  k  (long  wavelength)  <ok  pre¬ 
dicts  phonon-type  spectrum  whereas  for  k  large  (short 


of  free  particle.  These  results 
agree  qualitatively  with  the 
spectrum  first  proposed  by 


Landau3  and  verified  experi- 
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mentally  by  Henshaw  and  Woods  .  The  detailed  structure  of 
the  spectrum  (see  diagram)  for  intermediate  values  of  k 
depends  on  the  particular  form  of  the  interaction. 
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CHAPTER  7.  LIFE-TIME  OF  A  QUASI-PARTICLE 


The  interaction  between  the  excited  states  of  a 
many-boson  system  and  the  ground  state  gives  rise  to  a 
finite  life-time  for  the  excited  quasi-particles.  Con¬ 
sidering  the  ground  state  as  a  source  for  these  excita¬ 
tions  we  can  describe  the  coupling  between  the  source 
b^|0>  and  the  excitation  cT|o>  by  the  Hamiltonian 

H  =  l  03j,b^bk  +  ftc*^c  +  l  W(k)  [b£c  +  c"^bk]  •  (7.1) 

k  ~  ~  ~  k 

*v  <v 

where  W(k)  represent  the  quasi-particle-quasi-particle 
interaction  and  ft  is  the  energy  of  the  excited  state 

4. 

c  |0>.  The  equations  of  motion  for  b^  and  c  are 


ih>k  =  w^.b^.  +  W(k)c  ,  (7.2) 

ic  =  ftc  +  l  W(k)bk  .  (7.3) 

k  ~ 

Let  x(t)  and  ^T_(t)  denote  the  following  amplitudes 

X (t)  =  <0 | c (t) |1>  ,  (7.4) 

and 

ipR (t)  =  <0|bk  (t)  1 1>  .  (7.5) 

Taking  the  expectations  values  of  Eqs.  (7.2)  and  (7.3)  we 


find 


■ 
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iik(t)  =  03kipk(t)  +  W(k)X(t)  ,  (7.6) 

i-X(t)  =  ^x(t)  +  l  W  (k)  ip,  (t)  .  (7.7) 

k  ~  ? 

We  solve  these  equations  by  the  following  method:  Let 


*k(t)  =  I  PU)<M?)  e  i(S+!2)t  , 

C 

x(t)  =  l  P  (c )  x  ( C )  e‘i(?+S2)t  , 

C 

where  P(£)  is,  as  we  shall  show  below, 

p(c)  =  x(c)  (t=o)  +  J  ^k (s)^k(t=o) 

k 

Substituting  Eqs.  (7.8)  and  (7.9)  in  Eqs. 
we  find  a  set  of  algebraic  equations 


(7.8) 


(7.9) 


(7.10) 


(7.6)  and  (7.7) 


(c  -  <sk)^k (C)  =  W(k)x(C) 


(7.11) 


CX (C)  =  I  W(k)^k (c) 
k  ~ 


where 


(7.12) 


(7.13) 


We  can  think  of  Eqs.  (7.11)  and  (7.12)  as  eigenvalue 


equations.  In  matrix  notation  we  can  rewrite  these 
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equations  as 


[c-<$k  o  0  .....  -W  (k  i  ) 

M 

0  ®  .  -W(k2) 

ib 

k2 

-W(ki)  -W  (k2 )  .  Z 

N  -V  ^  J 

Non-trivial  solutions  exist  if  and  only  if 


0 


(7.14) 


W2(k) 

det  A  =  n  (C  -  6.  )  (C  -  l  - — )  =  0  (7.15) 

k  ~  k  c  -6, 

or  if 


W2  (k) 

MC)  =  e  -  I  - —  =  o  .  (7.16) 

*  ^'6k 

Equation  (7.16)  determines  the  eigenvalues  £.  The  eigen¬ 
functions  x(S)  an(3  ^(C)  are  real  (since  6^,  W(k)  and 

<v 

hence  e  are  real) ,  orthogonal  and  may  be  assumed  normalized: 


X  (C)  X  <C  '  )  +  I  ^(0^(0  =  6(t(.  (7.17) 

I  X2(C)  =  1  ;  I  1>q(C)1>p(C)  =  6q<p  •  (7.18) 

e  e  ~  ~ 

From  Eqs .  (7.8)  and  (7.9)  we  have 


' 
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i|>k  (t=0) 
X (t=0) 

From  these  two 


=  l  P(?)ipk(C) 
c  ~ 

=  l  P(c)x(c) 

c 

equations  one 


can  derive  Eq . 


(7.19) 

(7.20) 

(7.10)  . 


Now  using  Eqs.  (7.17),  (7.11),  (7.12)  and  (7.16) 

one  obtains 


x  W(k) 

x  (c )  =  1  <K.(c)  =  - ~  ,  (7.21) 

/A'(C)  ~  ( C  —  <5  k )  /A  '  (C  ) 


where 

•  .  w2  (k ) 

A  '  (C)  =  -A~  =  1  +  l  - ^  >0  .  (7.22) 

k  (C  "  6k)2 


If  we  now  assume  that  our  initial  conditions  are 


X  (t=0 )  =  1  ;  bk(t=0)  =  0 


(7.23) 


then 


P(C)  =  x(c) 


(7.24) 


and  we  have 


ipk  (t) 


W (k)  e  1 (c+fi) 
(C~6k)  A  '  (C)“ 


(7.25) 


and 
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x(t) 


e-i  (c+fi) t 

~TrTo 


It  is  convenient  to  write  x ( t )  as 


x(t) 


-iftt  -izt 

—5 — -  ( p  dz  ; — r 

2tti  r  A  (z ) 

a 


(7.26) 


(7.27) 


where  the  contour  a  contains  all  roots  £  of  A(z)  =  0. 
Clearly  the  residues  at  the  poles  z  =  £  of  the  integrand 
of  (7.27)  are 


2  it  i 


e-iCt 
A'  (S) 


We  shall  now  assume  that  the  eigenvalues  5  are  dense  on 
the  x-axis  and  evaluate  (7.27)  along  the  contour  a  shown. 
Hence 


V. 

s 

> 

s 

x(t)  = 


iftt 
e _ 

2iTi 


/  dx 


-ixt 


[ 


1 


A  (x) 


A+  (x) 


(7.28) 


where 


A+  (x) 


,  •  \  ...  ,  ,  .  0,  r  k  W  (k ) 

lim  A(x±in)  =  lim  (x±in-— 2  f 

^->0  n-^0  2  it  0  k 


dk ) .  (7.29) 
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Note  that  we  have  transform  the  sura  in  (7.16)  to  an 
integral  over  k.  Thus  we  obtain 


A+ (x)  =  x  -  A (x)  ±  ir  (x) 


(7.30) 


where 


and 


A  (x)  = 


n 


2tt' 


oo  2  9 

P  j  "JUIT.  dk 


o 


X 


(7.31) 


r  (x) 


JL 

2  TT 


P  /  k2W2  (k)  6  (x-6,  )dk 
0  K 


(7.32) 


Since 


1 _ 1_ 

A_(x)  A+(x) 


2ir  (x) _ 

(x-A (x) ) 2+  r2  (x) 


(7.33) 


we  have 


x(t) 


-ifit 

e 


TT 


dx 


-ixt 


r  (x) 


+  r 


(7.34) 


As  a  first  approximation  in  evaluating  (7.34)  assume  A 
and  r  are  independent  of  x.  Then  the  integrand  has  poles 
at  x  =  A±iF.  Choosing  a  contour  in  the  lower  half  plane 
we  get 


x(t) 


-i  (fi+A ) t  -rt 
e  e 


(7.35) 


Thus  A  represents  the  shift  in  the  energy  of  x  because  of 
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the  coupling,  and  T  is  the  life-time  of  this  state. 

These  quantities  are  obtained  from  the  coupled  equations 

Re  A (A  -  iT)  =  0  ,  (7.36) 

Im  A (A  -  i r )  =  0  .  (7.37) 

Therefore  A  and  r  can  be  determined  if  W(k)  and  <5,  are 

.K 

known . 


A  more  realistic  and  yet  solvable  Hamiltonian  for 
this  system  would  be  the  following 

H  =  ftc+c  +  l  wkb^bk  +  l  W(k) [bkc  +  cfbk] 
k  ~  ~  ~  k  ~  ~ 

+  |  I  U(k,j)btb  .  (7.38) 

k,j  ~  - 

In  this  Hamiltonian  one  takes  into  account  the  quasi¬ 
particle  self  interaction. 

In  the  next  two  chapters  we  shall  show  how  the 
two  quasi-particle  potential  W(k)  can  be  calculated  from 
v(k-j)  using  the  S-matrix  approach  and  the  G-matrix 


method . 


' 
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CHAPTER  8.  INTERACTION  BETWEEN  TWO 

QUASI-PARTICLES 


In  the  previous  chapters  we  treated  the  many- 
body  problem  via  the  equation  of  motion  approach. 

Here  we  shall  try  to  treat  the  Hamiltonian  (3.1) 
directly . 

We  define  the  creation  and  annihilation  operators 

» 

for  quasi-particles  by  3^  and  3k ,  with  the  condition 
that 


r  ft 

LKk' 


ft+i 


5. 

k/P 


(8.1) 


To  lowest  order,  we  write  the  operator  ak  as  a  linear 

combination  of  3^.  and  3jt  for  all  values  of  k  including 

"13 

the  zero-momentum  state 


ak  =  h'  (NBk  +  g'  (N6-k 

ag.  =  h 1  (k)  +  g'  opek 


(8.2) 


Here  h'  (k)  and  g'  (k)  are  real  functions  of  k  and 


/ 


h'  (k)  =  h'(-k) 


g'  (k)  =  g’(-k)  . 


(8.3) 
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For  the  transformation  (8.2)  to  be  canonical  we  have 


h 1  2  (k)  -  g  1  2  (k)  =  1 


(8.4) 


We  determine  h'  and  g'  by  substituting  (8.2)  in  the 
Hamiltonian,  rearranging  all  the  terms  in  the  normal 
order  and  equating  the  coefficient  of  the  term 
+  3k8_k  to  zero.  Thus 

^  ~  'v 


[ek-  y  +  A'  (k)  ]h'  (k)g'(k)  +  -  A  '  (k)  [h  '  2  (k )  +  g  '  2  (k )  ]  =0  , 

(8.5) 


where 


A'(k)  =  i  l  [  v  (0 )  +  v(k-£)]g'^(£)  =  A 1  (— k )  ,  (8.6) 


£ 


and 


A'  (k)  =  i  v (k-£ ) h '  (£) g 1  (£)  =  A’  (-k) 


(8.7) 


£ 


Using  Eqs .  (8.4)  and  (8.5)  it  can  be  shown  that  h'  and 

g'  satisfy  the  same  equations  as  h  and  g  which  we  defined 
in  previous  chapters  (in  the  lowest  order) .  Thus  we  shall 
drop  the  primes  on  h  and  g. 


The  Hamiltonian  written  in  terms  of  the  quasi¬ 
particle  operators  is 

H  =  e  +  l  Ek3k6k  +  Hx  ,  (8.8) 

k  ~  ~  ~ 


■ 
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where 


=  I 

k 


{  Uk  -  U  +  j  A(k)]g2(k)  +  |  A(k)h(k)g(k) 


(8.9) 


and 


Ek  =  [ek  "  y  +  A(k)Hh2(k)  +  g 2  (k )  ]  +  2A  (k)h(k)g(k) 


=  {  [£k  -  y  +  A  (k )  ] 2  -  A2(k)}^  . 


(8.10) 


The  last  relation  is  obtained  by  using  Eqs.  (8.4)  and 
(8.5).  The  chemical  potential  y  is  to  be  determined 
from  the  equation 


n  =  I  Nk  =  l  <0|a+ak|0>  =  £  g2(k)  .  (8.11) 

k  ~  k  ~  ~  k 

The  residual  interaction  expressed  in  terms  of  the 

3^ ' s  contains  all  the  quartic  terms  of  creation  and 

+ 

annihilation  operators  3^  and  3 y>  i.e. 


HI  =  ^  .1  6j«,m+n  [Fl^B-£em6n  + 

j£mn  ~  ~  ~  ~  ~  ~  ~  ~  ~  ~  ~ 


n 


+  F3Bje£B-nem  +  F4BK6-mB-n  +  F5B- j- '  <8‘12) 


where 


F1  =  2v(j-m)g  (£)h  (n)  [h  (j)h(m)  +  g(j)g(m)] 


(8.13) 


■ 


, 
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F2  ~  2  fv  ( +v  ( j+^)  ]  h  ( j  )  g  (£)  h  (m)  g  (n) 

+  v  [h  ( j  )  h  (£)  h  (m)  h  (n)  +  g  ( j  )  g  ( i )  g  (m)  g  (n)  ]  (8.14) 

F3  =  2v  ( j-™)  h  (^)  g  (n)  [h  ( j  )  h  (m)  +  g(j)g(m)]  (8.15) 

F4  =  v ( j-m)h ( j )h (£)g (m)g (n)  (8.16) 


F5  =  v  ( j-m)  g  ( j  )  g  (£)  h  (ra)  h  (n) 


(8.17) 


The  equation  of  motion  for  the  quasi-particle  operator 
•j* 

3^  can  be  found  from  the  Hamiltonian  ( 8.8 ), 

6g(t)  =  i  Ejlt)  +  i  I  (q,B  (t)  )  (8.18) 

where 


I  (q,3  (t) )  = 


-  y 

^  j£m 


q+ j , £+m 


v(q-f)  :  [h(q)a^a^a^+g(q)a^a_£am] 


(8.19) 


In  Eq.  (8.19)  the  symbol  :  implies  that  after  the  a's  are 
replaced  by  the  3's  from  (8.2),  then  the  resulting  expres¬ 
sion  must  be  rearranged  in  the  normal  order  with  respect 
to  the  3's. 


Now  consider  the  scattering  of  two  quasi-particles. 
Let  p  and  q  be  the  initial  momenta  and  m  and  n  the  final 
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momenta  of  the  two  quasi-particles.  The  scattering 
matrix  S  for  this  process^  is 


Smn;qp  ~  <?u}  out  I  9P  in  > 


p  tin  i 

<mn  out  3  p  out> 

1  q  1  ~ 


4* 

lim  <mn  out  I  3T (t) 


p  out  >  e 


-iE  t 

q 


(8.20) 


2  9 

Using  the  well-known  LSZ  formalism  ,  we  find 


(S-±; 


00 

f 


mn ;  qp 


,  r  .+  +  _iEat 

■  j  at<mn| (3  (t)-iE  3  ( t ) ]  jp  out >  e  c 

— oo  *3  q  q 


00  -iE  t 

=  -i  J  dt<mn  out |  I  (q ,  3  (t )  )  |p  out  >  e  ^ 


(8.21) 

where  we  used  (8.18)  to  get  the  last  relation.  The 
lowest  order  contribution  to  the  scattering  matrix  comes 
from  the  term  I(q,  3°Ut(t))  in  (8.21).  This  functional 
is  obtained  by  iterating  Eq.  (8.18)  once,  i.e.,  since 


hout  =  i  e  6+out  , 

q  q  q 


(8.22) 


therefore 


ef  ~  i  E  8  3+  +  il(q,  3OUt(t)  ) 
Mq  q  q  ~ 


(8.23) 


57 


This  is  a  valid  approximation  when  the  interaction  is 
weak.  With  this  approximation  we  can  calculate  the 
matrix  elements  in  (8.21)  and  do  the  time  integration 
with  the  result  that 


(S-l) 


mn  ,  qp 


-27riS  (E  +E  -E  -E  )V(pq;mn) 
m  n  <1  P  ~  ~  ~  ~ 


(8.24) 


where 

=  ^6P+q,m+n^v  +g  (S)g  (?)  1 

x  [h  (p)h (n)+g (p)g (n) ] +v (q-n)  [h (q)h (n) +g (q) g  (n) ] 


x  [h  (p)  h  (m)  +g  (p)  g  (m)  ]  +v  (q+p)  [g  (p)  h  (q)  (g  (m)  h  (n) 

+  h (m)g (n) ) +h  (p) g  (q)  (h (m)g (n) +g (m)h (n) ) ] }  .  (8.25) 


Thus  the  V-matrix  is  just  the  Born  term  for  elastic 
scattering  of  two  quasi-particles.  Had  we  used  the 
exact  equation  of  motion  in  (8.21),  the  V-matrix  would 
have  been  replaced  by  the  T-matrix. 


In 


the  barycentric  coordinate  system 


30 


P  =  k  =  -q 


m  =  k '  =  -n 


(8.26) 


and  the  V-matrix  becomes  (Eq.  (8.25) 
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V(k\k)  =  — ^  {  [v(k-k!  )+v(k+k' )  ]  [h(k)h(k')+g(k)g(k')]2 
~  ~  (2tt)  ~  ~  ~  ~  ~  ~ 

+  4v(0)h(k)h(k')g(k)g(k')}  ,  (8.27) 

where  we  went  over  from  box  normalization  to  the  continuous 
case.  Doing  partial  wave  analysis  we  find  that 

V  (k'k)  =  — {v  (k',k)  [l+(-l)£]  [h(k)h(k')+g(k)g(k')]2 

Xy  /  r.  \  D  Xy 

(2u ) 

+  4v£  (0  )  h  (k)  h  (k  '  )  g  (k)  g  (k  1  )  }  .  (8.28) 

For  | k |  =  | k ' |  we  have 

2 

V  (k)  =  — (v?  (k)  [1+  (-1)  Z]  +  A~4"  [v£(0)  +  (l+(-l)  £)v£(k)]}. 

*  (2ti)  Ek 

(8.29) 

Using  Eqs.  (3.30),  (5.6)  and  (5.12)  we  have 

=  4h2  (k)  g2  (k)  =  4  [  1  +  ^  F(k)]  |  F(k) 

Ek 

=  4  [1  +  pF  (k)  ]  p  F  (k)  ,  (8.30) 

where  we  defined  the  density  p  =  Thus  in  the  high 

energy  limit  if  we  assume  that  v£(k)  goes  to  zero  faster 
than  F (k)  then  we  have 

V„(k)  r — *  — “  r  p  F(k)v,(0)  .  (8.31) 

£  k^“  (2tt)  3 


, 


1  for  k 
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In  the  low  energy  limit  assuming  that  F(k)  = 
small  we  have 


A 


k->0 


if  p  >>  1 
if  p  <<  1 


(8.32) 


and 


Vk) 


k->0 


- v  [4p2(l  +  (-1)  ^ )  ]  v  (0)  for  p  >>  1 

(2tt)J  * 


— r  [5(1  +  (-1)^)  +  4  ]  v  0  ( 0 )  for  P  <<  1 

(2tt  )  36 

(8.33) 


Once  we  know  the  V-matrix,  (which  approximates 
the  T-matrix)  we  can  find  the  effective  potential  W(r) 
between  two  quasi-particles  from  the  equation 


T (q 1 ,q) 


iq-r  o 

=  jy/e  ~  ~  W(r)  ijj  (q  '  ,r)  d  r  , 


(8.34) 


T(q',q)  can  also  be  written  in  terms  of  the  phase  shifts 
(k)  as 


2m 

2 


T(q'  / q ) 


"ft 


2iq 


I 

£ 


(2£+l) P  ( 


q-q 


2i6 


£  qq 


■)  (e 


£ 


-  1) 


(8.35) 


Thus  in  our  approximation  Eq.  (8.35)  relates  the  phase 
shift  in  each  partial  wave  to  V.  Knowing  6^(k)  for  all 
values  of  k  enables  us  to  obtain  an  expression  for  the 


- 


two  quasi-particle  interaction  W(r),  according  to  the 
well  known  method  of  the  inverse  scattering  problem. 
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CHAPTER  9 .  THE  G-MATRIX 


In  this  chapter  we  develop  a  different  approach 
to  the  many-body  problem.  Suppose  that  the  Hamiltonian 
(3.1)  after  diagonalization  can  be  written  as 


H  =  E  +  l  wka£ak  .  (9.1) 

k  ~  ~  1 


In  this  equation  we  have  neglected  higher  order  terms 
V  ft 

such  as  l  s . .a. a. a- a.  etc...  Again  the  quasi-particles 
ij  1  3  1  3 

are  characterized  by  their  energy-momentum  relation  03k . 

We  divide  H  into  two  parts,  Hq  and  H  ^ , 


H  =  H  +  Hx 
o  I 


(9.2) 


with 


+  I 


EA&k 


(9.3) 


The  quantities  e,  Ek  and  H  are  given  by  Eqs.  (8.9), 
(8.10)  and  (8.12)  respectively.  Since  E  and  e  are  c- 
numbers  they  can  be  absorbed  in  H  and  Hq  respectively. 

We  would  like  to  point  out  here  that  the  specific  split¬ 
ting  of  H  we  chose  is  not  essential  to  this  method  and 
any  other  choice  of  Hq  and  HI  is  acceptable. 


For  a  given  state,  with  momentum  k,  the  difference 

to  -  E  measures  the  energy  shift  caused  by  the  interaction, 
k  k 


' 
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and  since  E^.  is  known  we  will  calculate  in  terms  of  E 

To  this  end  we  expand  the  operator  a  a  in  terms  of  M  ; 

<2  P  ?  P 


aa  =33  +  j  K(p,q;u,v)33  +  .... 

SB  5?  U,v£p,q  B-  ?! 


(9.4) 


Both  H  and  Hq  are  diagonal,  the  former  in  a^.  1  s  and  the 
latter  in  3^’s.  Using  (9.1)  and  (9.3)  we  find  the 
following  commutation  relation: 


[H,  a  a  ]  =  -(w  +  id  )  a  a 

p  q  p  q  p  q 


(9.5) 


[H  .  3  3  ]  =  -  (E  +  E)33 

O  PUPV  U  V  U  V 


(9.6) 


With  the  help  of  expansion  (9.4)  we  have 


[HI'  “qV  =  [H~H°'  “qV 


=  (Eq+VVWp+  £  K  <P  'S  ! "  'Y >  (Eu+Ev"VV  6u*V 

~  ~  ~  ~  ~  u,v^p,q 

(9.7) 


The  kernel  K  can  be  found  by  taking  the  matrix  element  of 
Eq.  (9.7)  between  the  state <0 |  and  bJVj0>,  where  x  and 
y  are  assumed  to  be  different  from  p  and  q. 


-1 


K(p,q;x,y)  (Ex+Ey_Wq_a)p)  <0  I  [HI  'apaq^x^y  I  0>  *  (9,8) 


i  X 
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Substituting  (9.8)  into  (9.4)  we  have 


a  a  =3B  +  — 

q  P  qPp  2 


l 


<0  I  [H  ,a  a  ] e  +  B+ I 0> 
I  q  p  u  v1 


u  /  v^p  ,  q  E  +  E  -  co  -  co 
~  ~  ~  ~  u  v  q  p 


;uev  •  0.9) 


We  now  define  the  V  and  the  G-matrix  in  terms  of  their 
matrix  elements 

<mn|v|pq>  E  -<  0  |  [H^. ,  eqBpJ  |  0>  (9.10) 

and 


<mn  G  pq>  =  -<0 


[H  ,a  a  ]  3^3^* 
1  q  p  m  n 


0> 


(9.11) 


Substituting  (9.9)  into  (9.11)  we  have 


<mn  G  pq>  = 


^  <mn  I  V I  u,  vxu ,  v  |  G  I  p  ,q> 

<mn|v|pq>  +  j  J 

u , v^p , q 


co  +  co  —  E  —  E 

q  p  u  v 


(9.12) 


This  integral  equation  relates  the  G-matrix  to  the  effec¬ 
tive  interaction  between  quasi-particles  V  (W  in  Chapter 
7).  We  note  that  (9.12)  differs  from  the  T-matrix  integral 
equation  in  two  ways.  First  instead  of  the  usual  two  boson 
potential,  we  have  an  effective  potential  which  depends  on 
the  properties  of  the  system,  such  as  the  number  of  exci¬ 
tations  and  their  distribution  etc .  The  second 

difference  is  in  the  energy  denominator  of  (9.12).  Here 


■ 


the  energies  of  intermediate  states  have  different  depen¬ 
dence  on  the  momenta  than  the  energies  of  the  initial  or 
final  quasi-particles. 
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Using  the  explicit  form  of  (Eq.  (8.12))  in 
(9.10),  we  find  that  the  matrix  elements  of  V  are  the 
same  as  given  in  (8.25).  Therefore  if  the  effective 
potential  is  weak  then 


<mn  G  pq>  ~  <mn|V|pq> 


(9.13) 


and  the  matrix  elements  of  S-l  are  proportional  to  the 

mnf  v-i  v  ol  nman  +-  o  r\- F  O. 

AUUi  -L-  _i_  ^  N —  ' —  *11' —  N— 


The  energies  co  and  oj  in  (9.12)  are  related  to 

q  P 

the  diagonal  elements  of  G  and  can  be  determined  in  the 
following  way:  First  we  calculate  the  commutator 


[«pVhJ  =  2“pap“p 


=  [a  a  , 

p  p 


+  HI] 


(9.14) 


Then  we  find  the  matrix  elements  of  (9.14)  between  the 

JL 

state  < 0 1  and  3p$pl0>  i*e., 


V0|Vp6X10>  =  <0 


[a  a  , 

P  P 


y  e  $,]  o+3+ 1  o> 

£  k  k  k  p  p1 


(Hi'aPyepBp|0> 


<0 


(9.15) 


- 
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Using  Eq.  (9.9)  we  find 


<0  la  a  3  +  3f I 0>  =  2 
P  P  P  P1 


(9.16) 


and 


<0I  i“pV  l  =  4E 


,t„  + 


(9.17) 


By  substituting  these  results  in  (9.15)  we  obtain  the 

following  relation  between  oo  and  E 

P  P 


wp  =  Ep  +  4<ppIgIpp>  • 


(9.18) 


T  Ti  rro 
—  **■  w 


neral  the  elements  of  G — matrix  depend,  on  ths 


momentum  of  the  center  of  mass  of  the  two  interacting 
quasi-particles.  If  we  assume,  however,  that  the  motion 
of  the  center  of  mass  gives  a  small  contribution  to  the 
G-matrix,  we  can  reduce  the  integral  equation  (9.12)  and 
then  correct  the  result  for  the  case  where  the  center  of 
mass  is  not  at  rest.  In  the  barrycentric  coordinates 
system  we  have  p  =  -q  =  k  and  m  =  n  =  k',  therefore  (9.12) 
can  be  written  as 


<k '  | G I k>  =  <k  '  I  V  I k>  + 


I  I 

P/k 


<k  '  |  V  |  pxp  |  G  |  k> 
"k  ■  EP 


(9.19) 


where  we  have  assumed  that  .  It  is  convenient  at 


<  -  .  *r 
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this  point  to  go  over  from  a  discrete  set  of  states  to  a 
continuous  set  by  defining 

G(k',k)  =  — - — rr  <k  '  I  G  I  k  >  (9.20) 

~  ~  (2tt)  ^  ~ 


V(k'  ,k)  =  — -  ■■  <k'  | v I k>  . 

~  ~  (2tt)J  ~ 


We  have 


,  V(k '  ,p)G  (p,k) 

G(k',k)  =  V  (k '  ,k)  +  i/dp  — - — ~  ~  (9.21) 

~  ~  '  ~  4  “k  "  Ep 

*V  *V 

where  V(k',k)  is  given  by  Eq.  (8.27).  We  observe  that 
the  denominator  of  the  integrand  in  (9.22)  may  vanish  for 
certain  values  of  p,  say  pQ,  which  are  different  from  k. 
Therefore,  we  assume  that  can  have  a  small  imaginary 
part.  Since  defined  by  (8.10)  is  a  real  quantity  the 
integral  in  (9.22)  will  be  finite  but  G(k',k)  will  also  be 
a  complex  quantity.  Let  us  write 

wk  =  \  "  2  Yk 

G  (k ' ,k)  =  5  (k' ,k)  -  i  T  (k* ,k)  . 

Substituting  for  o>k  and  G(k',k)  in  Eq.  (9.22)  and  separat- 
ing  the  real  and  the  imaginary  part,  we  find  two  coupled 
equations 


(9.23) 

(9.24) 
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o  V(k' ,p)  [H  (p,k)  (fik-E  )  +  jr  (p,k)yk] 
E(k',k)  =  V(k',k)+  |/d3p  — —  -  -  ~  ~  5  4  ~  ~  ~ 


(fik  -  v2  +  yk/4 


(9.25) 


r  (k '  ,k)  =  jf d3p 


V (k '  ,p)  [r  (p,k)  (fik-E  ) -yk  H(p,k)] 


<  °k  -  EP)2  +  ^/4 


(9.26) 


Also  from  Eq.  (9.18)  it  follows  that 


nP  =  EP  +  tt-  2(P'?> 


(9.27) 


and 


(9.28) 


Defining  F'(k',k)  by 


F'  (k' ,k)  =  r  (k‘ ,k)/r  (k,k) 


(9.29) 


and  substituting  into  (9.26)  ,  we  obtain 


(2tt  ) 


F' (k‘ ,k) =  x/d3p 


V  (k  '  ,  p )  [F  '  (p,k)  (^k“Ep) - 4^ -  5  ] 


-  EP)2  +  ^k/4 


(9.30) 


For  k ' =  k. 


F'  (k,k)  =  1 


and  therefore  (9.30)  becomes 
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V(k,p)  [F'  (p,k)  (SJk-E  )  -~p—  E(p,k)] 
t/<33P  - 5-^ - 5 -  =  1  .  (9.31) 

-  y  +  tk4 

Thus  we  have  three  unknowns  H(k,k'),  F'(k.k')  and  y^, 
that  can,  in  principle,  be  determined  from  the  three 
equations  (9.25),  (9.30)  and  (9.31).  Once  these  quantities 

are  calculated  the  energy  shift  5!^  and  the  life-time  of  the 
quasi-particle  y^  can  be  found. 
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CHAPTER  10.  CONCLUSIONS 

We  have  tried,  in  this  work,  to  develop  different 
approximation  techniques  for  the  treatment  of  a  many— boson 
system. 


As  we  discussed  in  the  introduction  and  in  Chapter 
2,  the  conventional  perturbation  methods  should  not  be 
used  to  treat  such  a  system.  This  is  because  the  S-matrix 
and  probably  the  transformations  that  diagonalize  H  are 
non-analytic  in  the  coupling  constant.  Thus  the  need  for 
a  systematic,  self  consistent  approximation  theory  remains. 

The  non-perturbative  theory  of  Bogoliubov,  when 
proposed,  seemed  promising  because  it  predicted  qualita¬ 
tively  the  correct  observed  results.  However,  there 
remain  many  difficulties  with  the  theory  such  as  the 
origin  of  the  singularities  of  the  transformations  employed, 
the  treatment  of  the  condensate,  high  order  corrections 

etc .  In  Chapters  3  and  4,  we  showed  how  one  can 

generalize  or  "derive"  a  generalized  Bogoliubov  transfor¬ 
mation  from  first  principles.  While  one  can  write  such 
an  expansion  consistently  the  resulting  equations  (aside 
from  lowest  order)  are  very  complicated  indeed.  The 
lowest  order  equations  obtained  from  the  generalized 
expansion  can  be  solved  consistently  (as  was  done  in 


' 
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Chapters  5  and  6)  without  separating  a  priori  the  zero 
state  in  the  Hamiltonian.  The  results  obtained  are 
essentially  the  same  obtained  by  Bogoliubov. 

A  more  practical  approach  to  the  problem  is  the 
model  Hamiltonian  approach.  In  Chapter  7,  we  showed 
how  to  construct  exactly  solvable  model  Hamiltonians 
which  are  relatively  simpler  to  deal  with.  In  order 
that  the  model  Hamiltonian  shall  describe  approximately 
the  many-boson  system  we  have  to  determine  the  effective 
interaction  between  two  quasi-particles.  While  this  in 
itself  is  a  difficult  problem,  we  have  shown  in  Chapter 
R  how  one  can  find  (approximately)  this  effective  potential 
via  the  S-matrix  method. 

Finally  in  Chapter  9  we  developed  yet  another, 
totally  different  approach  based  upon  a  'pairing'  type 
of  boson  operator  expansion.  This  expansion  yields  a 
complicated  integral  equation  for  our  G-matrix.  We  have 
shown  the  connection  between  the  G-matrix  and  the  effec¬ 
tive  interaction  of  quasi-particle  and  indicated  how  the 
quasi— particle  life-time  and  energy  shift  can  be  calculated 
using  the  G-matrix. 


'■  . 

- 
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APPENDIX  A 


We  want  to  investigate  the  S-matrix  as  a  function 
of  N.  To  this  end,  consider  the  Hamiltonian 


Using  the  same  procedure  as  of  Chapter  2  we  have 


(n),o  <hn)  S'  ^Jo> 


N 


N 


7T  n  c 

•  -j  n  . 
3=1  3 


°°  -  ( l  +  4ig )  q  . 

/  e  3  Hn.  (3-i)d3i 

■  OO  j  -J  -J 


(A. 2) 


Here  >  =  ip  is  the  unperturbed  wave 

^ny  n^  ,  ^2  r  *  *  *  * 

function  of  (A.l)  and 


n_,  -\ 


C  =  [/rT  n  .  !  2  J] 
n .  i 

u  J 


(A. 3) 


We  note  that  for  n^  =  odd  integer,  the  integrand  in  (A. 2) 

is  an  odd  function  of  q .  and  S ,  »  vanishes  for  this 

(n)  ,o 

case.  Hence  n.  must  be  even.  Let 

3 


n .  =  2m . 
3  3 


m . 
3 


0,1,2,... 


(A. 4) 


Then 
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(n)  ,o 


=  TT 


* : 
j  =  l 


'2m  . 
3 


«>  -Z2q2 

I  e  3  H2m(qj)dqj 

■  on  *•' 


N  /(2m.)!  _ 

~  n  - A -  (Z  -  1) 

„N  .  _  m. 

Z  3=1  2  3(m J! 


m . 
3 


(A. 5) 


with  Z  =  (1  +  4ig)^  .  The  probability  of  transition 
between  the  state  0  and  (n)  is 


p  ,  .  =  Is  .  .  I  2  =  Is  I2  n 

(n) r o  1  (n),o‘  1  0,0* 


N  (2m . ) !  m . 


(ft  D 


j  =  l  (m. !) 


2  4 


(A. 6) 


Here 


is 


0,0 


_N 

_M , 9  9  9 

z  “r  =  (i  +  i6g~)  “ 


(A. 7) 


and 


y  =  16g2 [ 1  +  16g2 ]  1 


(A. 8) 


Consider 


J  P  ,  . 

(n)  (n)'° 


.  N  oo  (2m .  )  !  m . 
2  v  J_  (Z)  3 


=  |  s  I  n  l  - 

0,0  j=l  m  ^  =  0  (irm ! ) 


2  4 


(A. 9) 


=  S 


o ,  o 


N 

n  (l  -  y)  2 

j=i 


s  |  (l  -  y) 

o  ,o  1 

N 


N 

‘2 


N 


(1  +  16g2 )  2  (1  +  16g2)2 


=  1 


(A. 10) 


(A. 11) 


(A. 12) 
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Hence,  while  the  dependence  on  N  of  each  q  might 

be  very  complicated,  the  sum  of  all  the  elements  Is ,  > 

1  (n)  ,o 

is  equal  to  1  and  is  independent  of  N. 


For  this  particular  model  we  can  calculate  the 
S-matrix  in  closed  form; 


S  (n)  ,  (m) 


S 

0,0 


N 

n 

j=i 


m .  +n  . 
2  3  D 


r iig_ 

^ l+4ig 


m .  +n  . 

_J _ 1 

2 


r  ( 


m  .  +n . +1 
3  3 


2 


) 


x 


2F1 


1-m . -n 
_ ] _ 

2 


1. 


l+4ig . 
8ig  ' 


(A. 13) 


with  n.+m.  =  even. 

n  n 

.j 
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APPENDIX  B 


Consider  the  Hamiltonian 


H  = 


N  to .  p  N  o 

l  ~Y  (c^i  +  ?i}  +  49  F(  l  qi)  6  (t) 
i=l  i=l 


(B .  1 ) 


We  want  to  find  under  what  conditions  the  S-matrix,  for 
the  system  (B.l),  will  be  non-analytic .  Thus  consider 


S  =  TT 

o ,  o 


N 

‘2 


/  aq- 


2  2  2  2 
dqN  exp [- (qx+ . .qN)  -  4igF  (q1+. .qN) ] 


(B .  2 ) 


Going  over  to  polar  coordinates  (see  Chapter  2)  we  have 
(N  =  k+2 ) 


0,0 


r(|) 


/ 

0 


-r^-4igF(r^)  k+1  , 

e  ^  r  dr 


r  (f) 


r  (|) 


2  -x  -4igF  (x) 
x  e  e 


dx 


0 


I 

n=0 


(~4ig) 

ni 


n  00 


I  dx  e 
0 


k 

-x  2 


x^  [F (x) ] n  dx.  (B. 3) 


Suppose  we  can  write 


F(x)  =  l  ax1  , 


£=0 


i 


(B .  4 ) 
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i  •  O  ^ 

then,  it  can  be  shown  that 


00  oo 

Fn(x)  =  [  l  a^xh n  =  l  C0(n)xZ 


£=0 


£=0 


£ 


(B .  5 ) 


Here 


n  1  m 

;  C_(n)  =  — —  V  (kn-m+kJa.C  . 

L  k  m-k 


o 


m 


ma 


o  k=l 


(B .  6 ) 


Hence  eq.  (B.3)  becomes 


1 


o°  /  .  n 

0,0  “  r  (N/2)  nJ=0  CILin)  r<|  +  £) 


(B.7) 


Define 


£=U 


B(n)  E  l  c„  (n)  F  (|  +  £)  . 


(B.8) 


Then 


S  (a)  =  _i _ y  (-4i-) 
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(B .  9 ) 


Using  the  ratio  test,  the  radius  of  convergence,  of 
the  series  (B.9)  is  given  by 


lim 

n->°° 


4  B (n+1 ) 
n+1  B  (n ) 


(B. 10) 


Thus  if 


lim 

n->°° 


(n+1) 
(n) 


(B.ll) 


then  the  S-matrix  is  a  non-analytic  function  of  g. 


